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ABSTRACT 

The Crab pulsar is well-known for its anomalous giant radio pulse emission. Past studies have concentrated only on the very bright 
pulses or were insensitive to the faint end of the giant pulse luminosity distribution. With our new instrumentation offering a large 
bandwidth and high time resolution combined with the narrow radio beam of the Westerbork Synthesis Radio Telescope (WSRT), we 
seek to probe the weak giant pulse emission regime. The WSRT was used in a phased array mode, resolving a large fraction of the 
Crab nebula. The resulting pulsar signal was recorded using the PuMa II pulsar backend and then coherently dedispersed and searched 
for giant pulse emission. After careful flux calibration, the data were analysed to study the giant pulse properties. The analysis includes 
the distributions of the measured pulse widths, intensities, energies, and scattering times. The weak giant pulses are shown to form a 
separate part of the intensity distribution. The large number of giant pulses detected were used to analyse scattering and scintillation 
in giant pulses. We report for the first time the detection of giant pulse emission at both the main- and interpulse phases within a single 
rotation period. The rate of detection is consistent with the appearance of pulses at either pulse phase as being independent. These 
pulse pairs were used to examine the scintillation timescales within a single pulse period. 

Key words, pulsars - neutron stars - emission mechanism - giant pulses 
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1. Introduction 

: >: 

, C*") Identified as the supernova remnant that resulted from SN 1054, 

■ ■ the Crab nebula is one of the strongest radio sources in the 
OO sky, and it harbours the young neutron star PSR B0531+21. 

The pulsar is visible across the entire observable electromag- 
' t^J" netic spectrum, and at radio wavelengths it is the second bright- 
\ ' est pulsar in the northern sky. P SR B0531+21 was discovered 

■ , by IStaelin & Reifenstein (1968), soon after the discovery of 
' ^ — 1 ■ pulsars. This pulsar is noted for several features including the 

I near orthogonal alignment of the magnetic and rotational axis 

■ • i-H that gives rise to the observed interpulse emission. The aver- 
\ ■ age emission profile of the pulsar, obtained by averaging the ra- 
, H \ dio emission from many rotations of the star, exhibits a number 
■- . . i of features that change qu ite remarkably with radio frequency 

dMoffett & Ha nkins 1994). The single pulses show a large vari- 
ation in amplitude and duration as a function of time. The most 
enigmatic of t hese are its occassional intense bursts kno wn as 
giant pulses dHeiles et alJ Il970t IStaelin & Suttonl Il970h . The 
giant pulses can be ex tremely narrow, of the order of 0.4 ns 
(Hankins & Eilek 2007) and the pulse flux can be several 1000 
times the average pulse flux. The ultrashort durations of the gi- 
ant pulses imply ver y high equivalent brightness temperatures 
dHankins et al.l [2003b indicating that they originate from non- 
thermal, coherent emission processes. In this work, we define 
giant pulses as the pulses with a significantly narrower width 
than the average emission and contain a flux of at least 10 times 
the mean flux density of the pulsar. 



The Crab pulsar is one of just a handful of pulsars that have 
been shown to have giant pulse emission. Some other pulsars, 
like the young Vela pulsa r, also show narrow, bursty emission 
called giant micropulses dJohnston et alJ l200lb . The fluxes of 
these micropulses are within a factor of 3 times the average 
pulse flux. In the pulsars that show giant pulse emission, the 
pulse intensity and energy di stributions exhibit power-law statis- 
tics dArgvle & Gowerlll972l) . w hile the giant mi cropulses give 
rise to log-normal disnibutions (Cairn s et alJl200ll) . In contrast, 
the bulk of the pulsar population have pulse intensities and en- 
ergies that follow either a normal or an expo nential distribution 
dHesse & Wielebinskill974l lRitchings 19761) . This indicates that 
the giant pulses and micropulses may form a different emission 
population. 

The Crab giant pulses have been studied by different groups, 
yet the nature of the emission process remains elusive. In the 
very early studies at low sky frequencies, the data were af- 
flicted by dispersion sm earing and scattering dHeiles et aljl970t 
Idow er & Argviel l 19721) . but the power-law nature of the inten- 
sity disu'i bution of giant pulses was identified. In the next ma- 
jor study. ILundgren et alJ d 1995b discuss a multi- wavelength ob- 
servation of giant pulse emission, and note the possibility of 
a weak giant pulse emission population at radio wavelengths, 
which they are unabl e to resolve owing to insufficient sensitivity. 
Sall men et alJ d 19991) found that the Crab giant pulses are broad 
band at radio wavelengths. They also determine giant pulse spec- 
tral indices in the range of -2.2 to -4.9 using their widely spaced 
observation bands and 29 simultaneously detected giant pulses. 
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Observations by Hankin set al.l (1200 3) revealed that giant pulses 
at 5.5 GHz contain nanosecond wide subpulses and the presence 
of such narrow features ha s been predicted in numerical mod- 
elling by Weatheralll dl998l) . At these frequencies the radio emis- 
sion character of the Crab pulsar changes, with the interpulse 
emission becoming dominant. A multi-wavelength radio obser- 
vation of Crab giant pulses with widely spa ced frequency band s 
(0.43 GHz and 8.8 GHz) is presented by ICordes etall (12004). 
who discu ss the effects of scintillati on ov er a wide range of fre- 
quencies. IPopov & Stappersl (120071) and lEilek et all d2002l) in- 
vestigated pulse wi dth distributions and find that narrow pulses 
tend to be brighter. lBhat et all d2008l) carried out a similar analy- 
sis in addition to scattering and dispersion variations in the neb- 
ula. All of these studies point to the peculiarity of the Crab pulsar 
and its puzzling emission process, and motivates further study in 
finer detail using a large number of pulses. For the work dis- 
cussed in this paper, we utilised t he wide band capabilities o f the 
new pulsar machine, PuMa-II ( Karuppusa my et al.l 120081) and 
the Westerbork Synthesis Radio Telescope (WSRT) in the coher- 
ent tied-array mode. At small hour angles, the synthesised beam 
of the WSRT effectively resolves out the Crab nebula, reduc- 
ing the nebular contribution to the system temperature. Thus the 
WSRT and PuMa-II combination makes this study much more 
sensitive in terms of signal-to-noise ratio achieved, and in num- 
ber of pulses than was possible in the past. The rest of the paper 
is organised as follows: in §2 we describe the observational set 
up and data reduction, flux calibration is discussed in §3, the 
giant pulse characteristics are discussed in §4. We report detec- 
tions of double giant pulses in §5, and the scattering analysis is 
presented in ^6] 

2. Observations and data reduction 

The radio observations of the Crab pulsar reported here were 
carried out as part of a multi-wavelength observation with the 
Integral y-ray telescope and the WSRT on 1 1 October 2005. The 
WSRT observations were from UTC 03>?56™5O! to 09 1 !36™20 s . 
with a break of three minutes in the middle of the observation 
to switch data disks. The results of the y-ray observations will 
be reported elsewhere. 

The pulsar was observed at eight different sky frequencies 
in the L-Band, which is the most sensitive front-end receiver at 
the WSRT (Tsy S = 30 K). The sky frequencies (see Table. [TJ 
were chosen to be free of radio frequency interference. Two or- 
thogonal polarisations of 8x20 MHz analogue signals from each 
telescope were 2-bit sampled at the Nyquist rate of 40 MHz. The 



Parameter Value 

Observation duration 21420 s 

Start Epoch 53654.726505 (MJD) 

Sky frequencies 1311 a , 1330, 1350, 1370, 1392" 

1410, 1428 ab , 1450 MHz 

Bandwidth 8 x 20 MHz 

Nominal T sys 30 K 

Beam size 21" x 1741" c 



a These frequencies are not uniformly spaced to avoid interference. 
b This band was not recorded due to disk failure. 
c The beam size varies as a function of the observation time. See text 
for details. 

Table 1. Telescope parameters and observation details. 
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Fig. 1. Total intensity of a coherently dedispersed giant pulse at 
the main pulse phase detected in all recorded bands at 4.1 /is 
resolution. The total dispersion delay of 24.9 ms across the seven 
bands was removed for this plot. The lower most panel shows the 
pulse after combining the signal in all seven bands. The pulses 
displayed here are scaled relative to the pulse at 1330 MHz. 

telescope was operated in the tied-array mode in which coherent 
sums of the sampled voltages were formed in dedicated adder 
units resulting in 6-bit summed voltages. A coherent sum was 
achieved by determining the instrumental phase offsets between 
the telescopes using observations of a strong calibrator source. 
These phase offsets, combined with the geometrical phase off- 
sets required for tracking the source are applied to each tele- 
scope. The resulting values were then read off as 8-bit data and 
recorded in the PuMa-II storage nodes. This resulted in a total 
of 13.5 Terabytes of raw data. After the observation, the data 
were processed offline using the open-source pulsar data pro- 
cessing software package DSPSPQ- A 32-channel synthetic co- 

1 http://dspsr.sourceforge.net/ 



R. Karuppusamy et al.: Giant pulses from the Crab pulsar 



3 




0.2 0.4 0.6 

Pulse Phase 



Fig. 2. The plot shows the average pulse profile (top panel) and 
the total intensity for six of the seven recorded bands in greyscale 
(lower panel). The striped nature of channels at 1330 MHz and 
1390 MHz comes from the overlap in the adajcent frequency 
bands. The roll-off of the filters used in the system is also seen 
as a reduced intensity at the band edges. A low-level extended 
feature is seen at the edge (also visible in the top panel as the 
elevated baseline in the right side of the main pulse) of each band 
which is due to the 2-bit quantisation noise and is only visible in 
long exposures. 



herent filterbank was formed across each 20 MHz band with co- 
herent dedispersion applied across each of the channels using 
the dispersion measure (DM) of the pulsar. We obtained the DM 
(= 56.742) from the Crab p ulsar ephemeris maintained by the 
Jodrell Bank ObservatorjQ ( Lvne et a"Dl 19931) at the epoch clos- 
est to our observation. Frequency resolution was preserved so 
that studies of spectral indices, scintillation, and scattering could 
be carried out. 

The total intensity was computed for each pulse from the 
dedispersed data. Giant pulses were detected by computing the 
peak signal-to-noise ratio (denoted by S /N). The giant pulse de- 
tection threshold was set at S /N > 7<x in each band, where <x is 
the off-pulse root-mean-square noise fluctuation. Pulses below 
the detection threshold were discarded to ease storage require- 
ments. The original sampling time was 25 ns. The 32-channel 
filterbank and the choice of 4. 1 yus final time resolution resulted 
in 8192 phase bins. The time resolution of 4.1 jjs was chosen 
to match the estimat ed scattering timescale available at the time 
( Sallmen et al. 1999). However, it is known from recent work by 
Bhat et alJ (120081) that single pulses at these radio frequencies 
can be as narrow as 0.5 /vs. In addition to the single pulses, aver- 
age pulse profiles with 128 frequency channels in each 20 MHz 
band were formed every 10 seconds. 

The reduced data consisted of ~2 1 000 giant pulse candidates 
in each recorded band. An example candidate is shown in Fig.Q] 
where the pulse was detected in all bands. In the offline analysis 
stage, these candidates were combined in software using only 
pulses that show the expected dispersion delay. This method 

2 http://www.jb.man.ac.uk/ pulsar/crab.html 



ensures that spurious signals were filtered out in our analysis. 
After combining in software, 12959 giant pulses were identi- 
fed to have occurred simultaneously at all observed sky frequen- 
cies. Of the 12959 pulses, 1 1384 were detected at the main pulse 
phase and 1370 at the interpulse phase of the average pulse pro- 
file. 

The data were folded and the single pulses were formed us- 
ing the DSPSR soft ware package and a poly nomial determined 
by using TEMPO dTavlor & Weisberplll989h . The folded pro- 
files formed in each 20-MHz band were combined in software 
to validate the DM used. The combined data are shown in Fig. [2] 
as a frequency-phase image and shows no smearing, confirming 
that the value of DM is correct. A similar procedure was used to 
combine simultaneous giant pulses in all seven bands. Some arti- 
facts of the 2-bit systems of the individual telescopes are visible 
once the profile is summed for the entire six-hour long observa- 
tion. The width of these artifacts match the dispersion smearing 
in the bands as seen in the top panel of Fig. [2] The quantisation 
noise is 12 % for a single telescope, whose signal is sampled us- 
ing 2-bits dCooper| [l970). Since signals from the 14 telescopes 
of the array were coherently summed, the uncorrelated quantisa- 
tion noise was reduced by a factor of Vl4. The resulting noise 
of 3.7% is considered too small to be problematic in the analysis 
that follows. In many stages of the analysis, extensive use of the 
PSRCHIVE dHotan et alJl2004 utilities was made to view and 
validate the pulsar data and to compute the S /N used in later 
analysis. 

3. Flux calibration 

To establish a flux scale for the observed giant pulses, the mean 
system flux needs to be computed. The mean system flux is pro- 
portional to the r.m.s noise variations at the te lescope outp ut and 
can be expressed by the radiometer equation (Dicke 1946), 

5* min — ) j ( Fj 

y/N p ■ B ■ T~ t 

where, S,„i„ is the r.m.s system noise in Jy, S svs the total sys- 
tem noise, N p the number of polarisations (= 2), B the band- 
width in MHz (= 140), and Tj„, the integration time in seconds. 
The total system noise in flux density units in eqn. (Q]i is the to- 
tal equivalent system temperature divided by the telescope gain 
(S sys = T tota i I G). For the WSRT, when signals from the four- 
teen 25-m parabolic dishes are combined, the resulting telescope 
gairQis G = 1.2 K.Jy -1 . The term T tota i, can again be expressed 
as 

Ttotai = T S y S + fy(t) ■ T cn . (2) 

The term T cn is the contribution of the Crab nebula to the system 
temperature, w hile f v (t) is a time-depen dent factor explained be- 
low. Following Biet enholz et al.l d 1997b . we express the total flux 
of the Crab nebula at frequency v (in GHz) as Scn = 955 v~ - 27 
Jy, from which T c „ is computed. The WSRT is an east-west ar- 
ray and the coherent addition of the telescope signals results 
in a 21" x 1741" fan beam. The Crab nebula is an extended 
source of size Q„, = 6' x 4', so the WSRT's fan beam resolves 
the Crab nebula in the east-west direction. This in turn reduces 
the nebular contribution to the T sys . However, the width of the 



3 The telescope gain is 1.34 KJy~' for an ideal array combiner. The 
reduction in gain is attributed to losses in the formation of the tied-array 
signal. 
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Fig. 3. The upper panel shows the change in minimum detectable 
signal S „„■„ within a 4. 1 //s time interval during the first 5 hours of 
the 6-hour observation. The hour angle of the source is displayed 
on the top ordinate axis. The lower panel is the observed change 
in peak signal-to-noise ratio of the detected signal. The depen- 
dence of the signal-to-noise ratio on the hour angle of source is 
discussed in the text. 

WSRT's fan beam is not a constant, but is a function of the ob- 
servation time. While the source is being tracked, the effective 
width of the synthesised beam changes with hour angle (HA) 
and it is expressed as Qa(0 = Oc« ■ A/D ■ cos(HA). In this ex- 
pression HA =t- RA, where t is the local sideral time, the max- 
imum baseline D = 2700m, and RA is the right ascension of the 
Crab pulsar. The fraction of the nebular contribution can be ex- 
pressed as f v (t) = Q,A(t)/Q, cn , which reaches its minimum value 
of 0.13 at zenith. As the source is tracked towards the horizon, 
the projected distance between the dishes decreases and £2^(f) in- 
creases. Consequently, the observing system becomes less sensi- 
tive toward larger hour angles, or when the source rises and sets. 
This time dependence of the system noise is included in our flux 
calibration. The variation in S m! „ is shown for a bandwidth of 
140 MHz, N p =2 and r = 4.1 fis in the upper panel of Fig. [3] A 
plot of the pulse intensity during the observation (lower panel of 
Fig-EJ confirms this reduction in sensitivity. 

The peak flux of the giant p ulses were computed using the 
modified radiometer equation (Lorimer, D. R. and Kramer. M. 
120051) for the pulsar case, S pea k — (S/N) ■ S„,i n . With the above 
considerations of the nebular contribution to T tota i and with 
T sys = 30K in the WSRT's L-Band, the system retained suffi- 
ciently high sensitivity in the first 15000 seconds of the obser- 
vation. Two other factors have been neglected in this calibration 
procedure and do not contribute significantly to the T sys : the rel- 



ative change in the orientation of the WSRT's fan beam and the 
Crab nebula over the course of observation and the partial shad- 
owing of three telescopes out of the 14 for HA > 54° (the last 3 
hours of our observation). 



4. Single-pulse statistics 

For the analysis that follows, all pulses that were flux-calibrated 
as described in the previous section were used. The discussed 
change in system sensitivity does not limit this analysis thanks 
to our careful flux calibration procedure. While approximately 
70% of the pulses were detected in all seven bands simultane- 
ously, the rest were detected in two or more of the seven bands 
recorded. For the results described below, where applicable, only 
those pulses that were detected in all seven bands were used and 
explicitly mentioned. 

4. 1 . Pulse intensity distributions 

The giant pulse fluxes of the Crab pulsar contribute to 
the long exponential tail of the single pulse intensity his- 
tograms ( Arg yle & Gowerl 1 1 972b . while the normal pulsars 
show Gaussian or expone ntial pulse intensity distributions 
dHess e & WielebinskJ Il974h . Fig. [4] shows the average pulse 
flux distribution for pulses detected in at least two of the seven 
recorded bands. The average pulse flux is computed by integrat- 
ing all emission within the equivalent width, W eq of the giant 
pulse (see §14.4I >. This value is averaged over the pulse period to 
obtain the average pulse flux. The pulse in each band was de- 
tected based on a threshold of Icr. A pulse detected in two bands 
satisfies the V2 X 7 = 9.89<x limit. In the first three hours of 
the observation (when the system was most sensitive), the flux 
equivalent system noise in 4.1 /is is 109 Jy. Averaged over the 
pulse period, a pulse of S/N = 9.89<x corresponds to an aver- 
age pulse flux density of 3.9 Jy. This implies that it is sensitive 
to all pulses greater than 27 x (F), where (F) = 14raJy is the 
average flux density of the Crab pulsar. Therefore, the flux dis- 
tribution computed here contains a good fraction of weak giant 
pulses compared to those reported elsewhere (see Table. [2]). 

The intensity distributions displayed in Fig.|4]shows at least 
two components: a peak at or below ~ 4 Jy - the weak pulses that 
may comprise the trailing part of the normal pulse distribution. 
The next component peaking at ~ 20 Jy resembles a lognormal 
distribution with a power-law tail. The bright giant pulses result 
in the extended power-law tail and is described by Nf °c F a , 
where Nf is the number of pulses detected in 1 .8 Jy flux intervals 
of F. The value of a = -2.79 + 0.01 and a = -3.06 + 0.06 was 
determined from the best fits to the data in the interval 1 1 8 Jy 
< F < 2000 Jy and 40 Jy < F < 596 Jy for the giant pulses in 
the main- and interpulse, respectively. Visual inspection of Fig. 
|4]shows that the distribution is multi-modal, with giant pulses in 
the region F ^ 10 Jy and the pulses below this limit possibly 
representing normal pulses. 

It is worth noting the differences in the intensity distribu- 
tions displayed in Fig. [4] While the distribution of the giant 
pulses in the main pulse phase shows a clear turn over at ~20 
Jy, the emergence of a bimodality in the region containing weak 
pulses is evident in the intensity distribution of the interpulse gi- 
ants. The distribution corresponding to the interpulse phase also 
shows a flattening in the 10-30 Jy region. The clear excess of 
weak pulses in both the distributions in the region F < 4Jy is 
due to our method of setting W eq = 4.1 yus (equal to the time 
resolution). In this case the emission window we considered is 
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Fig. 4. Distribution of the pulse intensity of all giant pulses 
detected at the main- and interpulse phases in the upper and 
lower panels, respectively. The long tail results from the giant 
pulse emission. The best fit power-law curve is shown with slope 
-2.79 + 0.01 for the pulses in main pulse phase and -3.06 + 0.06 
for the pulses in the interpulse phase. Both distributions show an 
excess near 4Jy and come from the rounding off in W eq . [see text 
for details]. 



Reference 


Frequency 


Threshold 




(MHz) 


(Jy) 


Lunderen et al. (1995) 


800 


120.0 


Ipopov & StappersTHoOTh 


1197 


5.9 a 


iBhat et all d2008f) 


1300/1470 


22.3 b 


This paper 


1373 


3.9 



a Equivalent average pulse computed flux from the quoted 6<x peak 
flux density of 142 Jy, assuming 0.036 pulse duty cycle. 

b Average pulse flux density extrapolated for 7cr threshold, 4.1//S 
time resolution and pulse duty cycle as 0.036. 

Table 2. Reported sensitivity to the Crab giant pulse observa- 
tions in the literature. 



dominated by noise or weak and narrow pulses. The slopes of 
the power-law models obtained here can be compared t o the val- 
ues reported earlier. Fig. 4 of Lundgren et al. ( 1995j) shows a 
slope of -3.46 + 0.04 for data at 800 MHz, which is slightly 
steeper than the slopes of t he ma in- and interpulse distributions 
derived h ere.lCordes et al.l (|2004l) derive a value of ~ -2.3 at 433 
MHz and Bha t et alJ (120081) found -2.33 + 0.14 at 1300 MHz, 
which are comparable to the slope the main pulse intensity dis- 
tribution in our work. The slopes of the intensity distribution re- 
ported here generally agree considering the effect of low number 
statisics and/or dispersion smearing in the observations reported 
elsewhere. While this experiment was sensitive to much lower 



fluxes, the long observation time has also enabled the detection 
of rarer bright pulses. 

4.2. Pulse energy distributions 
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Fig. 5. The cumulative probablity distribution of the energy in 
giant pulses detected at the main pulse and the interpulse phases 
in the upper and lower panels, respectively. The y-axis is the 
fraction of the total number of pulses and pulse energy is plotted 
on the x-axis. Also shown are the occurrence rates per minute, 
second and hour. 

The relative occurrence rates of giant pulses is displayed as 
a cumulative probablity distribution of the individual pulse ener- 
gies in Fig. [5] The pulse energy is computed by multiplying the 
equivalent width, W eq , and the average pulse flux. As described 
in 34.11 we computed the best fits to the cumulative probablity 
distributions of the main- and interpulse giants. The power-law 
curve with a = -2.13 ± 0.007 and a = -1.97 ± 0.006 fits the 
data for pulse energies at the main- and inter pulse phases, re- 
spectively. The break s een at ~2000 Jy.^s is con sistent with the 
break value reported bv lPopov & S tappers! ((2007). The emission 
at the interpulse phas e shows a somewhat shallo wer power-law. 

It is known from iPopov & S tappers! J2007) that the power- 
law index has a width dependence, varying from -1.7 to -3.2 
as the pulse width increases. Bas ed on this variation, the i ndex 
we find is in good agreement with Popov & Stap persl (120071) and 
IBhat et all (120081) (-1.88 ± 0.02 at 1300 MHz). However, we 
fit only a single power law unlike the two power-law fits found 
by these authors. Partial fits to the low-energy pulses yield more 
than two components, with shallower power-law indices indicat- 
ing a simple dual-component fit is insufficient. One explanation 
for this can be the bias introduced by setting W eq = 4.1 yus for 
narrow pulses, overestimating the pulse energy. However, this 
can only be a minor contribution and is an argument that there is 
a clear break in the intensity distribution. To compare the occur- 
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rence rates we see here, we proceed to derive the rates from the 
arrival times of the giant pulses in the next section. 

4.3. Giant pulse rates 

The distribution of the separation times between successive gi- 
ant pulses is plotted in Fig. [6] If the giant pulses are mutually 
exclusive events independent of each other, then the arriva l time 
separation follows a Poisson process (Lund gren et al.ll 995). The 
probablity of a giant pulse occurring in the interval x is then 
given by P(x) = /j.x.e^\ where /j. is the mean pulse rate. Since 
our data only consist of giant pulses, we expected to see an ex- 
ponential reduction in the separation time between the pulses. 
Fig.[6]shows the fits to the separation times at both the inter- and 
main-pulse phases. 




50 1 00 

Separation time (s) 

Fig. 6. The symbols show the distribution of separation times 
between successive giant pulses at the main- and interpulse 
phases and the solid lines are the best fits to the distribution. 
The top ordinate axis corresponds to the curve and data for the 
pulses at the main pulse phase and are offset by 450 for clarity. 

Functions with an exponential decay with time constants 
1/t = 1.1 + 0.02 and 1/t = 0.172 ± 0.003 are in excellent agree- 
ment with the data at the main- and interpulse phases, respec- 
tively. From the values of r, the mean giant pulse rates are one 
main- pulse giant every 0.9 seconds and one inter pulse giant 
every 5.81 seconds observed above our threshold limit of 3.9 Jy. 
At these frequencies, the interpulse giants are comparatively less 
numerous as is evident from our data. For comparision, the inter- 
pulse giants are brighter and more f requent in frequency bands 
above 5.5 GHz (ICordes et al.l [2004). The combined rate of the 
giant pulses (fit and data not shown) is one pulse every 0.803 
seconds. The foregoing discussion confirms earlier predictions 
that th e giant pulse rate increases with frequenc y for the Crab 
pulsar dLundgren et al.lll995tlSaiimen et al.ll l999h . The effect of 
the WSRT's sensitivity reduction towards the end of the obser- 
vation, as displayed in Fig. [3] may have contributed to the long 
tail of the distribution, where fewer pulses were detected than in 
the first half of the observation. However, the rate derived here 
is robust, since the system had sufficiently high sensitivity in the 
first half of the observation. 



4.4. Width distributions 

The equivalent pulse width, W eq is defined as the width of a top- 
hat pulse with height equal to the peak intensity of the pulse. W eq 
for the giant pulses detected in all seven bands was computed. 




pulse width, W gq (/2s) 

Fig. 7. Plot of intensity against pulse width for the main- 
and interpulse windows in the top left and lower left panels. 
Histograms of equivalent pulse widths are shown in the top right 
and lower right panels. The distribution has an exponential en- 
velope. For pulses with computed Weq < 4.1 /is due to random 
noise fluctuations the widths were rounded off to 4.1 [is. 

The results are displayed in panels on the right in Fig. [7] We 

express W eq as 



W. 



eq = j X V // X4.1 [IS, 

1 ma x 



(3) 



where I max is the peak intensity, /, the intensity in the pulse emis- 
sion window defined by bins i = n\ ■ ■ ri2 and is equal to 1ms in 
our case. Thus W eq can be viewed as the equivalent width of a 
rectangular pulse in [is that has the same area as the giant pulse, 
with height I max . 

The g iant pulses at these frequencies can be quite narrow. For 
instance, iBhat et al.l (120081) find pulse widths to be 0.5 [is and 
lEilek et all d2002l) found 0.2 /is. Our method of data reduction 
allowed a time resolution of 4.1 [is, so pulses with W eq < 4.1/is 
were taken to have a width equal to 4.1 [is. This results in some 
pulses being underestimated in flux and overestimated in equiv- 
alent width. The computed equivalent widths range from 4. 1 /is 
to ~120 [is, and we find that bright pulses tend to be narrow 
as seen i n the left hand pane ls of Fig. [7] This was also sug - 
gested bvlSallmen et al.l(ll999l) and shown bv lEilek et all d2002). 
IPopov & Stapp ers (2007) found a similar behaviour in addition 
to a width-dependent break in the power-law fits to the pulse- 
energy distribution. 

In the seven closely spaced radio bands observed, we note 
that a vast majority of the pulses have widths larger than 4.1 
[is. This is seen in the pulse width histograms at the two pulse 
phases, displayed in the panels on the right in Fig. [7] The dis- 
tribution shows a peak at ~16 [is, which is 4 times our ultimate 
time resolution in the main pulse, and the peak shifts towards 
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narrower timescales for the interpulses. We find less than 9% of 
the pulses with W eq =4.1 //s, indicating that the majority of the 
pulses show wider widths than our time resolution. The shape 
of the width distribution is similar at both the main- and inter- 
pulse phases. The contribution to the tail region of the distribu- 
tion comes from scatter broadened pulses. 

4.5. Spectral index of giant pulses 

The data were recorded in 7 different radio bands each 20 MHz 
wide in the frequency range 1300-1450 MHz, and several thou- 
sands of pulses were detected simultaneously in all bands. The 
spectral index of individual pulses was computed by modelling 
the flux variation of a giant pulse as S (v) oc v k . Here, S (v) is the 
flux of the giant pulse at frequency v, and k the spectral index. 
The histograms of the derived spectral indices are displayed in 
Fig. [8] for the giants at both pulse phases. A large dispersion in 
the spectral index is seen, with values -1.44 + 3.3 for the main- 
and -0.6 + 3.5 for the interpulse giants. 
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Fig. 8. Histogram of spectral indices for the giant pulses de- 
tected at the main pulse (bottom panel) and the interpulse phase 
(top panel). The spread in the distributions is indicative of fitting 
errors. See text for details. 

These spectral index values are quite a bit shallower than 
those detected previously (see Introduction) over wider fre- 
quency separations. We therefore consider the effects of diffrac- 
tive interstellar scintillation (DISS) on the spectral index esti- 
mates. Strong DISS results in pulse intensity variations within 
each of the seven bands. The effect of scintillation is to modu- 
late the observed pulsar signal in both time and frequency. This 
is seen as regions of enhanced or diminished brightness in a 
grey scale plot of the intensity as a function of time and fre- 
quency. These regions are known as scintles. We estimate the 
scintillation bandwidth based on the pulse scatter timescales, 
t s = 395 + 50/zs at sky frequency of 200 MHz, as reported in the 



work o f Bh at et all J2007). We further make use of their revised 
Tj oc v -3 ' 5 frequency scaling and consider that the scintillation 
bandwidth and scattering timescale are related by 2nAvjT s = C\ 
, where the consta nt C\ = 1.05 for a thin scattering screen 
(ICordes et al.l20 04). From these considerations Av,/ as 0.25-0.38 
MHz in the 1300-1460 MHz band. On examining a few giant 
pulses by eye, it was clear that some of the scintles are resolved, 
while some were narrower than our channel width of Av = 0.625 
MHz. Thus, in the flux obtained by integrating the signal in the 
20MHz-wide bands, the scintles tend to average out. This im- 
plies that scintillation does not cause the spread in the individual 
giant-pulse spectral indices. Moreover, with such narrow scin- 
tillation bandwidths, averaging over many giant pulse spectral 
index determinations as we have done here would give an aver- 
age spectral index that reflects the true average spectral index. 

Refractive interstellar scintillation (RISS) cannot corrugate 
the spectra of single pulses, since the pulse intensity variations 
due t o RISS are noticeable in observation of the order of a few 
days (Lundgren et aD ll995h . However, the pulses do have a sig- 
nificant structure that is intrinsic to the emission process. One 
example is displayed in Fig. \T\ and these pulses do contribute 
to the spread in the computed spectral indices. In this figure, it 
is clear that the leading short burst shows considerable varia- 
tion across the seven bands, while the scattered trailing part of 
the p ulse is correlated across frequency. This is again similar to 
what Hankins & Eilek (2007) find, as shown in their Fig. 4, but 
at a much higher frequen cy of ~9 GHz. 

ISallmen et al. (1999) find that the spectral index variation 
is between -4.9 and -2.2 based on 29 pulses they observed 
in two bands centred at 1.4 GHz and 0.6 GHz. The spread 
in the indices computed here and that of Sallm en et al.1 (l999) 
points to the stochastic nature of the giant pulse emission process 
and/or the disturbed plasma flow in the magnetosph ere caused 
by strong plasma turbulence (lHankins & Eilek 2007). The giant 
pulses used in this analysis were detected in all seven bands and 
represent 70% of all detected pulses in our data. Since each of 
our bands is 20 MHz wide, detection in seven bands implies an 
emission bandwidth of at least Av = 140 MHz. This suggests 
that the emission bandwidth of Crab giant pulses is potentially 
greater than Av/v = 0.1, unlike t he giant pulse emission from 
the millisecond pulsar B1937+21 (Pop ov & Stappersll2003l) . We 
note that the Av/v = 0.8 for the Crab giant pulses reported by 
Sall men et al.l(ll999l) was based on 29 simultaneous giant pulses 
from their 90-minute observation (~ 161086 stellar rotations). 
Those 29 pulses could have been chance detections, while the 
Av/v = 0.1 limit derived here comes from a much larger sample 
of giant pulses so is more robust. We detected a total of 17587 
giant pulses, of which approximately 4000 were detected in less 
than 7 bands. Clearly it is impossible to include the pulses de- 
tected in only a few bands in this analysis as that would increase 
the dispersion in the spectral indices computed; however, this 
lack of detection in all bands, for pulses which were clearly de- 
tected in the other bands, is an argument for there being some 
narrow band effects that appear to modulate the giant pulse in- 
tensity. 

5. Double giant pulses 

During direct inspection of some giant pulses, it was noticed that 
occasional giant pulse emission was evident at both the main- 
and interpulse phases within a single rotation period of the star. 
To determine how many such pulses were present, the follow- 
ing search algorithm was used. First, the giant pulses detected 
in all seven bands were combined in software across the fre- 
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quency bands. The pulses were then averaged over polarisation 
and frequency to create single pulse total intensity profiles. The 
search algorithm was made sensitive to emission at both emis- 
sion windows (main- and interpulse) by traversing each pulse 
profile twice; in the first pass, the emission peak and phase in- 
formation was recorded, following which a search is made in the 
other emission window i.e. if a pulse was detected at the main 
pulse phase we check whether a pulse is also seen at the inter- 
pulse phase. All pulses that show signal > 5<x in the second emis- 
sion window are collected separately. The pulses returned by the 
search procedure were examined by eye to validate the double 
pulse nature. To our knowledge, this is the first instance of this 
phenomena being reported. A total of 197 pulses that show emis- 
sion at both pulse phases were found in our data set above the 
5<x detection threshold. 
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Fig. 9. Detected double giant pulses shown as a ratio of the main 
pulse to the interpulse flux. The x-axis shows time since the start 
of the observation. 



To consider how likely this is to happen by chance, we note 
that the observation lasted 643263 rotations of the star and 11584 
and 1375 giant pulses were found at the main- and interpulse 
phases, respectively, above the Icr detection threshold in each 
band. Since these giant pulses were detected in all seven bands, 
the effective threshold is now V7x7o- = 18cr. If the 18cr crite- 
rion is used to search for the double pulses, a total of 17 pulses 
are seen. In other words, only 17 pulses in the 197 detected 
show S IN > 18cr in either of the two emission windows. Let 
the giant pulses occurring at the two pulse phases be indepen- 
dent events, with individual probablitites P(A) and P(B). The 
chance of two giant pulses occurring within a single rotation 
period is the joint probablity P(A,B) = P(A).P(B). Thus the 
chance of detecting a giant pulse above the 18cr threshold limit 
at the main- and interpulse phases are P(A) - 11584/643263 
and P(B) = 1375/642263 leading to P(A,B) = 3.5 x 1(T 5 . We 
therefore expect a total of P(A, B) x 643263 = 24 pulse periods 
with pulses at both phases in our data. The detection of 17 pulses 
is thus consistent with the expected 24 pulses. 

As seen above, combining the seven bands improves sen- 
sitivity and allows the detection of weaker pulses. Considering 
pulses with S /N greater than 5cr in the second emission win- 
dow resulted in the detection of an additional 180 double pulses. 
While the 197 pulses detected are not sufficient to perform mean- 
ingful statistics of these pulses, in § 16. II we use our population 
of double giant pulses to study scintillation and scattering within 
a 0.5 rotation of the pulsar. 

Although the appearance of the pulses in the same rota- 
tion period is consistent with the individual occurrence rates, 



we compared the GP properties at each phase. In the double 
pulses, the emission in the interpulse phase is typically narrower 
(W eq S= 16 [is) than the emission at the main pulse phase and 
pulses at the main pulse phase are typically brighter, as shown in 
Fig-E] In both cases this is consistent with the known population 
of GPs at each phase. A similar analysis to the one in {34.3l was 
done to determine the rate of double pulses and a rate of 1 pulse 
in 84 seconds, or one in 2545 rotations of the star was found to 
have giant pulse emission at both pulse phases. Thus, given the 
narrowness and very low occurrence rates of these pulses, they 
were easily missed in earlier observations. 



6. Single-pulse scattering 

The frequency resolution and large bandwidth of our data ben- 
efits scattering and scintillation checks on the individual pulses 
in two ways. First, the pulses detected in 7 bands are combined 
in software to give 224 channels across the 140 MHz bandwidth 
allowing examination of scintillation. Second, the large band- 
width of the combined pulse increases sensitivity and makes it 
possible to identify low-level extended scatter tails. To charac- 
terise the scattering time t s in the pulsar signal, we computed 
the extent of pulse broadening in the individual giant pulses. If 
the pulses are scattered by a thin-screen between the source and 
the observer, the pulses can then be modelled as a n one-sided 
expon ential with a vertical rise and a rapid decay (Williamson 
1972). This can be written as 



f(t) = 



e -t(T, 




if t > 
if t < 0. 



(4) 



This model was fit to the data using a least-squares min- 
imisation and the 1/e time derived from the models was taken 
as t s of an individual g iant pulse. It is known from the work 
of ISallmen et al. (1999), that a single one-sided exponential is 
not sufficient to model the complex structure of the giant pulses 
at this frequency. However, The large majority of pulses in our 
data show that the single exponential model agrees within 10% 
error. Therefore, we proceeded with the single exponential fits. 
The values of t s as a function of observing time and their dis- 
tribution are shown in the upper and lower panels of Fig. [10] 
respectively. The reduction in the scattering time towards the 
end of the observation is consistent with scattered pulses tending 
to be dimmer, hence below the detection threshold. Only suffi- 
ciently bright pulses are detected in the sensitivity limited part 
of the observation, as discussed in §[3] The scatter tail is also 
not discernible from the system noise in this part of the observa- 
tion, limiting the determination of t s . However, there were fewer 
pulses so they did not contribute to the distribution of t, (lower 
panel of Fig. ITOb significantly. 

The lower panel of Fig. [10] shows an exponential envelope 
in the distribution of t s . The individual pulse scattering time 
varies from 4.1 [is to ~90 [is. The large number of pulses in 
the distribution with r s « 4 [is is related to our ultimate time 
resolution of 4.1 [is.. This also implies that a large fraction of 
the pulses have scattering time t s < 4.1 /j . s. At a slightly ear- 
lier epoch than our observations, Bhat et al. (2007) determined a 
value of t s = 395 + 50 [is at 200 MHz. Using their revised fre- 
quency scaling of t s oc y -3-5±0.2 > tne sca ttering time at the centre 
of our band (1373 MHz ) is 0.47 ± 0.05 [is. At a slightly later 
epoch, iBhat et al.1 d2008l) find a value of t s = 0.8 ± 0.4 [is at 
1300 MHz, which contrasts with the value of 8;«s at 111 MHz 
(or 1.4 [is at 1300 MHz using a v~ 3 5 scaling law) reported by 
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Fig. 10. Upper panel: a plot of the values of time constant r s 
from the fits to the scattering in the entire duration of obser- 
vation for those pulses detected at the main pulse phase. The 
lower panel displays a histogram of the time constants obtained 
from exponential fits to the scatter tails of individual pulses. The 
values of r, < 4. 1 fj.s are from pulses narrower than our time 
resolution and those that are likely to be free of scattering. 



iKuzmin et al.l (120081) . With our data, we are not sensitive to scat- 
ter times below 4. 1 ps, but to the dispersion seen in the histogram 
of scatter times in Fig. [TOl shows that variations can even be ex- 
pected within a single observation of six hours. We again refer 
to Fig. Q] for an example of the extreme form of this variation: 
the different parts of the same pulse show different scattering ef- 
fects, imp arting a significant struct ure to the pulse. In their work 
on DISS. ICordes & Rickettl (jl998) emphasise that considering 
the 1/e time equal to t s is only valid for a thin screen and does 
not always hold. In light of the limited validity in interpreting 
the 1/e time and the spread in the values of scatter times found 
in our analysis, we suggest that the scattering in the direction 
of Crab pulsar cannot be modelled by single thin screen. The 
spread in t s ranges from < 4.1/zs to ~ 120/zs in our ~ 6 hour- 
observation. This proves most of the scattering cannot be due to 
the ISM, as the line of sight through the ISM does not change 
rapidly enough to explain these variations. Therefore, the bulk 
of scattering should orginate in the Crab nebula. The nebula can 
clearly give rise to a complex screen or changes in the struc- 
tures in the vicinity of the pulsar that give rise to the short-term 
changes in scatterin g time ( Ba cker et al.ll2000l iLvne et al.ll200H 
Sall men et al.l fl999). The scattering of pulses cannot be in the 
pulsar magnetosphere. In that case the pulses at lower frequen- 
cies that originate higher up in the magnetosphere should show 
lower scatter times, because according to the standard pulsar 
models, the number d ensity of charged particles i s lower in the 
upper magnetosphere (Lvubarskii & Petrov afl998l) . However, t s 



scales with frequency as v~ 3 - 5 (Popov et al. 2006), and this does 
not support the hypothesis that scattering could have its orgins 
in the pulsar magnetosphere. 

The diffractive scintil lation timesc a le, At nrss at this fre- 
quency was estimated by ICordes et al.l (120041) as 25.5s, based 
on pairs of single pulses with sufficient S /N. However, the pulse 
pairs they used were separated in time by a few pulse periods. 
Since our data has good frequency resolution (224 frequency 
channels across 140 MHz), and we detected several pulses with 
multiple components, we proceeded to estimate possible varia- 
tions in the scintillation time on shorter timescales. 



6.1. Scintillation within single pulses 

The scintillation timescale within single pulses was estimated 
using those pulses that show well separated components and the 
double pulses discussed in ^5] The search for at least two com- 
ponents in single pulses was carried out based on the component 
separation of ~25 /is. This was done by examining the pulses by 
eye, after an automated first pass. The first pass provided 45 1 gi- 
ant pulse candidates, 368 of those displayed at least two distinct 
shots in the main pulse phase, and 18 candidates were found in 
the interpulse phase. The 197 double pulses were included in 
this analysis. Assuming that the two shots of pulses are intrin- 
sic to the pulsar emission and that the scattering screen remains 
stable within a pulse period, any scintillation would affect the 
two components similarly, introducing a correlated frequency 
structure. The scintillation timescale is then the 1/e point along 
the time axis of the 2-dimensional intensity correl ation function , 
C(6v, t) = (I(t, v).I(t + r, v + 6v)} of the spectrum dCordesll 986). 
The computed correlation coefficients between the two compo- 
nents and the double pulses are displayed in Fig.fTTI 

The correlation coefficient of ~0.4 for many pulse compo- 
nent pairs is in exc ellent agreement with the value derived by 
ICordes et al.l (|2004). They derive a value of 0.33 considering the 
giant pulses to be 100% polarised, amplitude modulated, scintil- 
lated shot noise. It also implies that these components have un- 
dergone similar scintillation effects, ruling out the possibility of 
any variation in the scattering medium on these timescales. The 
average correlation coefficients computed for the double pulses 
is consistent with the average value computed for the widely 
spaced pulse components (pulses in the top panel of Fig. ITTb . 
Since a clear roll-off in the values of correlation coefficient is 
not seen in the data presented here, we conclude that the scintil- 
lation tim escales are longer th an 14 ms, which is entirely consis- 
tent with ICordes et all (120041) . 

7. Discussion 

To our knowledge this is the largest collection of high time- 
resolution giant pulse analysis presented in the literature. 
Even though some features of the giant pulse emission like 
the giant nano shots a re in the process of being explained 
(lHankins & Eilek 2007), several questions still remain about the 
pulsar emission mechanism in general and the giant pulse phe- 
nomena in particular. From the measured pulse widths and the 
observed structure in many pulses, it is evident from the analysis 
presented in this paper that the giant pulse emission is a manifes- 
tation of temporal plasma changes in the pulsar magnetosphere. 
The observed giant pulse rates are further evidence for this tem- 
poral variation, because if the mechanism responsible for the gi- 
ant pulses is active on timescales longer than a pulse period, a 
clear excess of giant pulses separated by a single rotation period 
can be expected. On the basis of the giant pulse arrival times, it 
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Fig. 11. Correlation coefficients of the spectra within a single 
pulse period. Top panel shows correlation between the two com- 
ponents of giant pulse, while lower panel is the double giants. 
The separation between the components t is shown in the ab- 
scissa. 



was concluded that the observed giant pulse emission does not 
come fr om a steady emission beam loosely bound to the stellar 
surface (Lund gren et al.lll995t[Sallmen et al.lll999l) . We confirm 
that our data do not support such a model, for if such a beam 
with random wobbles operates, a characteristic width in the gi- 
ant pulses can be expected. In other words, the distribution of the 
pulse widths would be normally distributed with a mean width. 

The power-law n ature of the gia n t puls e intensity distribu- 
tions was shown by Lundgren et al. (1995), and they inferred 
that the normal pulses formed a separate part of the intensity 
distributions. In this work, we have shown conclusively that 
the giant pulses consist of two distinct populations especially 
for those pulses found at the inter pulse phase. We see a def- 
inite change in the shape of the distribution of pulse energies 
as we go to lower energies and we also see a slight broadening 
of the pulses. These pulses still seem to be distinct from what 
might be called "normal pulses": they are still narrower than 
most subpulses and are at least 27 times brighter than the nor- 
mal pulses. The slope of the distribution containing these pulses 
is different from rest of the intensity distribution. These pulses 
could possibly be the t railing part of the distribution inferred by 
iLundgren et alj (119951) . Moreover, how these relate to the pre- 
cursor emission is unclear, which can clearly be improved upon 
using the double giant pulses. While there is evidence of a broad- 
ening of the pulses as they weaken in intensity, they do not ap- 
pear to be as broad as standard s ubpulses. This fi nding has im- 
plications in the model derived by Petrova (2004), where a clear 
power-law distribution is explained, but not a weak giant pop- 
ulation. The power-law index derived also has implications for 



interpretin g giant pulse e mission on the ba sis of self orga nised 
criticality dBaket al.ll 19871) . as suggested bv lCairnsi (120041) . 

The spectral index of the Crab giant pulses reported in this 
work suggests that the emission bandwidth is at least Av/ v > 0. 1 
and may approach the upper limit Av / v _ 0.2 predicted in nu- 
merical models by IWeatheralll (119981) . lHankins & Eilekl (120071) 
find a similar emission bandwidth at 9.5 GHz. Moreover, the 
average spectral index of giant pulses at the interpulse phase is 
flatter than the giant pulses at the main pulse phase. This possi- 
bly explains the dominant and bright nature of interpulse giants 
at v > 5 GHz. We note the prominent emergence of bimodality 
in the intensity distribution of the interpulses relat ive to the main 
phase pulses. Furthermore, (lHankins & Eilekll2007l) find upward 
drifting emission bands in the spectrum of the interpulses giants 
and not in the main pulse giants. These differences strongly sug- 
gest a different nature to the inte rpulses. To explain the drifting 
emission bands. ILvutikovl (120071) derived an excess plasma den- 
sity of ~10 5 and a large Lorentz factor of the emitting particles 
of the order of ~10 7 , and this condition is satisfied close to the 
light cylind er over the mag netic equator. However, the model 
proposed by ILvutikovl d2007l) is only valid for v > 5 GHz, where 
the emission bands are observed. While results from our obser- 
vations can neither support nor rule out this model, the differ- 
ence in pulse intensity distributions we find indicates that the 
interpulse giants are different in nature. 

It is worth noting that the pulsar signal is a stochastic process 
that contributes to the measurement noise of the pulsed intensity. 
This is especially true in the case of giant pulse emission, where 
pulsed flux can exceed 1500 Jy, an order of magnitude greater 
than the system equivalent flux density (SEFD) of approximately 
145 Jy. Source-intrinsic noise increases the measurement uncer- 
tainty of various derived parameters, such as the pulse d flux den- 
sity, p ulse width, scattering time, and spectral index Ivan Stratenl 
(2009). In addition, any temporal and/or spectral correlations - 
either intrinsic to the giant pulse emission or induced by inter- 
stellar scintillation - will also affect the uncertainties of any de- 
rived parameters. The vast majority of the pulses presented in 
this analysis have average flux densities that are lower than the 
SEFD, and we do not expect that self-noise will significantly al- 
ter the results of this analysis. To accurately quantify the impact 
of self-noise on parameter distributions (such as those presented 
in Figures 4,5,7, and 8) would require extensive simulations that 
are beyond the scope of the present work but may provide addi- 
tional insight in a future paper. 

The previously unreported double pulses we found are con- 
sistent with the occurrence rate on a purely probabilistic ba- 
sis. Collecting even more of these pulse pairs would allow for 
better checks of the statistics of occurrence to ascertain that 
they are chance occurrences and not indicative of some longer 
term underlying phenomenon driving the giant pulse emisision. 
Moreover detecting more of these pulses at higher time resolu- 
tion would provide furth er insight into the nature of these pulses. 
lHankins & Eilekl (|2007) found that the giant pulses at the inter- 
pulse phase show an additional dispersion when compared to 
the pulses at the main pulse phase. The closest pulse pair they 
were able to examine were separated by 12 minutes. One may 
gain new insight into the excess dispersion seen at the interpulse 
phase by examining the double giant pulses, which are the clos- 
est giant pulse pair possible. 

Scattering analysis of single pulses presented in this paper 
show a va riety of scattering t imes and corroborates with the 
analysis o flSallmen et al.l(ll999l) . They show that scattering from 
multiple screens or a single thick screen is excluded because of 
the observed frequency independence of the pulse component 
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separation. From this it was concluded that the multiple compo- 
nents that make up the giant pulses are intrinisic to the emission 
mechanism. Using multiple components and the double pulses, 
we conclude that the scintillation timescales are greater than 14 
ms, which indicates that there are no large changes in the number 
density of the scattering medium along the line of sight through 
the nebula on similar timescales. That the multiple components 
we detect in the giant pulses are spaced by at least 25 //s im- 
plies that the magnetosphere and/or the plasma does not change 
on these timescales, if the source intrinsic emission is less than 
25 //s. On the other hand, giant pulses may consist of overlap- 
ping nano shots. In this case the competing models make use of 
plasma turbulence leading to modulational instablity ( Weat heraiH 
119981) or the induced Compton scattering of low-frequency radio 
waves dPetroval2 004) in the magnetosphere to explain the origin 
of the nano shots. While with our data we are not sensitive to the 
pulses less than 4. 1 fi& duration, there is an indication that the 
emission bandwidth Av/v > 0.1, suggesting that the pulses can 
potentially have structure as narrow as 3.6 ns at this frequency. 

8. Conclusions 

The large collection of single pulses we gathered has allowed us 
to perform a range of statistics with the data. After careful flux 
calibration, a detailed analysis of the pulse intensities, energies, 
widths, and separation times was done by computing distribu- 
tions of these quantities. In the single-pulse intensity distribu- 
tions, we find a clear evidence of two distinct populations in the 
giant pulses. The giant pulse separation times show a Poission 
distribution, and the rate of occurrence of giant pulses was deter- 
mined. Spectral indices for a large number of giant pulses were 
computed with the narrowly spaced multi band data. Significant 
dispersion in the spectral indices was found and a small negative 
average spectral index was found for the main- and interpulse 
giants, and they are flatter than the average pulse emission. We 
also note that in some cases there is evidence for intensity modu- 
lation with bandwidths that are smaller than the full band but not 
consistent with scintillation effects. The previously undetected 
double giant pulses were presented and we find that they are 
not more frequent than would be expected by chance. The scat- 
ter time for a large number of giant pulses was determined by 
modelling the scatter broadening as an exponenial function and 
the distribution of scatter times was computed. The double giant 
pulses were reported for the first time and it is found that they 
are not very different from the normal giant pulses. Using multi- 
ple emission components either at the main- or interpulse phase 
and the double giant pulses, we find no evidence of variation of 
the scattering material on timescales shorter than 14 ms based 
on the correlation coefficient computed for emission within a 
single-pulse period. 
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